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Abstract 

Classical integrable Hamiltonian systems generated by elements of the Poisson 
commuting ring of spectral invariants on rational coadjoint orbits of the loop algebra 
(2,-R) are integrated by separation of variables in the Hamilton-Jacobi equation 
in hyperellipsoidal coordinates. The canonically quantized systems are then shown 
to also be completely integrable and separable within the same coordinates. Pairs 
of second class constraints defining reduced phase spaces are implemented in the 
quantized systems by choosing one constraint as an invariant, and interpreting the 
other as determining a quotient (i.e., by treating one as a first class constraint and 
the other as a gauge condition). Completely integrable, separable systems on spheres 
and ellipsoids result, but those on ellipsoids require a further modification of order 
0{'h?) in the commuting invariants in order to assure self-adjointness and to recover 
the Laplacian for the case of free motion. For each case - in the ambient space i?", 
the sphere and the ellipsoid - the Schrodinger equations are completely separated in 
hyperellipsoidal coordinates, giving equations of generalized Lame type. 



Introduction 



A general method for realizing integrable Hamiltonian systems as isospectral flows in 
rational coadjoint orbits of loop algebras was developed in [AHP, AHH1-AHH4]. This 
approach begins with a moment map embedding of certain Hamiltonian quotients of sym- 
plectic vector spaces into flnite dimensional Poisson subspaces of the dual fll(r)+* of the 
positive frequency part of the loop algebra gl(r) (or certain subalgebras thereof). The Adler- 
Kostant-Symes (AKS) theorem [A, K, S] then implies that the spectral invariants provide 
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commuting integrals inducing isospectral flows determined by matrix Lax equations. The 
level sets of these commuting invariants are shown to determine Lagrangian foliations on the 
rational coadjoint orbits, and hence completely integrable systems. Finally, a special set of 
canonical coordinates, the spectral Darboux coordinates are introduced, in which the Liouville 
generating function, which determines the linearizing canonical transformation, is expressed 
in completely separated form as an abelian integral on the associated invariant spectral curve. 
The resulting linearizing map is essentially the Abel map to the Jacobi variety of the spectral 
curve, thus providing a link, through purely Hamiltonian methods, with the algebro-geometric 
linearization methods of [Du, KN, AvM]. This approach has been applied to the study of 
a large number of integrable classical Hamiltonian systems, as well as the determination of 
finite dimensional quasi-periodic solutions of integrable systems of PDE's [H, HW, AHH3, 
AHH4, W, TW]. 

In the present work we focus on the case gl(2)"'"*, taking an equivalent approach to 
integrability based, first of all, on separation of variables in the Hamilton- Jacobi equation. 
The relevent "spectral Darboux coordinates" in this case simply reduce to hyperellipsoidal 
coordinates. The purpose of this reformulation is to prepare the passage to the corresponding 
quantum systems and the study of integrability and separation of variables in the associated 
Schrodinger equation. As it turns out, each such classical integrable system has an integrable 
quantum analogue, for which the Schrodinger equation is completely separable in the same 
coordinates. One case of separation of variables in such systems, - the quantized Neumann 
oscillator (an anisotropic harmonic oscillator constrained to the surface of a sphere) - was 
studied in [BT] , and the results extended to the quantum Rosochatius system in [Mc] . Other 
special cases, involving quantized free motion in various symmetric spaces and reductions 
thereof, were studied in [K, KM, KMW, BKWl, BKW2, ORW, Ku, To]. All these 
systems may be placed in a loop algebra setting using the moment map embedding of [AHP] , 
and canonically quantized. The resulting formulation is equivalent to a Gaudin spin chain [G, 
Ku], with 5u(2) replaced by 0^(2), and the separation of variables interpreted as a "functional 
Betheansatz" [Ski, Sk2]. 

In [Mo] , the algebraic geometry of a number of classical integrable systems constrained 
to quadrics in R"^ was examined. The integration of these and related systems was given a 
loop algebra formulation based on (2, C) and reductions thereof in [AHP, AHH4]. In 
the present work, such systems will be reexamined in terms of separation of variables in the 
Hamilton-Jacobi equation. Their quantum analogues will then be studied through constrained 
canonical quantization, making use of the loop algebra formulation to identify the commuting 
invariants in terms of "quantum determinants". The corresponding Schrodinger equations 
will be shown to separate within the same coordinates as the classical systems. Constraints 
leading to dynamics on spheres and ellipsoids will be shown to lead to integrable quantum 
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systems, also separable in the same coordinates as the classical ones. 

In Section 1, the appropriate loop algebra formulation of the systems in question is given. 
In each case, the coadjoint orbit is identified with the quotient R^"/(Z2)", and integrable 
isospectral flows are examined both in this space and on constrained submanifolds identified 
with the cotangent bundle of a sphere 5""^ C R"^, or an ellipsoid S'^~^ C R". The key step 
consists of using a Lagrange interpolation formula to express the invariant spectral polynomial 
in terms of its values at the associated spectral divisor points, and noting that these values 
coincide with the squared canonical momentum components. The separation of variables 
follows from identification of residues in the interpolation formula. 

In Section 2, the corresponding quantum systems are obtained by canonical quantization 
in the ambient phase space R^"^ before quotienting. The resulting Schrodinger equation is 
again expressed in terms of hyperellipsoidal coordinates through Lagrange interpolation, and 
the completely separated form is deduced, again by identification of residues, giving various 
types of generalized Lame equations. The associated one dimensional Schrodinger operators 
are seen to give the quantized form of the invariant spectral curves (cf. [Sk3]). In the case 
of the constrained systems on the sphere 5""^ the same formulation, together with second 
class constraints, leads without difficulty to completely separable integrable systems. In the 
case of the ellipsoid £'^~^, however, a new problem arises, since the separation of variables, 
while holding for the Schrodinger equation, does not hold for the volume element, giving 
rise at first to non self-adjoint operators. This is easily rectified by noting that the resulting 
Schrodinger operators are nevertheless self-adjoint with respect to a scalar product determined 
by a different measure than the one associated to the induced volume form. Conjugating the 
operators by the map relating the two scalar products, self-adjoint operators are obtained 
with respect to the standard measure. However, to recover the Laplacian in the case of free 
motion, a further scalar term of order 0{fi^) must be added to the quantum Hamiltonians. 
Since this correction gives the same semi-classical limit and does not destroy the integrability 
(or separability), it may be viewed as a satisfactory quantized version of the associated classical 
integrable systems. 

1. Classical Systems in 0l(2)+* 
la. Ambient Space. 

Following the general approach of [AHP, AHH2, H], we define a Poisson map 

j;,:R2n.^ ~[(2)+* (Lla) 
Xi:(x,y)^Aro(A) (Lib) 
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where 

Hi 




5^ XiVi-m v^n 



En 
- ..^ i=l A-Oi I , (1.2) 

= 1 A— Oi Z_(i=l A — Oi 

X, y G R"^ have components (a;^, 2/j)i=i,...„ and {/^i, Q;i}i=i,...n is a set of 2n arbitrary real 
constants, the CKj's being chosen as distinct. Here, the loop algebra 0l(2) consists of smooth 
maps X : — > from a fixed circle S^, centred at the origin of the complex A-plane, to 
0((2), the subalgebra 01(2)"'" consists of elements X e 0((2) admitting a holomorphic extension 
to the interior of S^, and the (smooth) dual space ^[(2)+* is identified with the subalgebra 
of elements X admitting a holomorphic extension outside S^, with X{oo) = 0. The 
dual pairing < , > is defined by integration: 

< X,Y > :^ (f ti(X(X)Y(X))dX (1.3) 
27rz J SI 

X e J'r^*(2) ~ 0r(2)-, Y e 0[^(2). 

To assure that A/q G 0[(2)+*, we must choose the constants {ai} entering in the definition (1.2) 
of the Poisson map J a to be in the interior of 5"^. The image of the map is a coadjoint orbit 
in 0l(2)+*, identified with the quotient R^"^/(Z2)'^ of the phase space R^" by the symplectic 
action of the group (Z)2 of refiections in the coordinate hyperplanes, and {xi,yi}i=i^,,,n are 
canonical coordinates defining the standard symplectic form 



Fixing an element 
we let 

XiX) :=y+jVo(A) 



uj = dxi A dyi. (1-4) 



Y={1 _^Je.Z(2,R), (1.5) 



2 Mi 



= {r - J + 9 1 ^ ^ 0^(2). (1.6) 

The ring X^^^ of commuting invariants is chosen, according to the Adler-Kostant-Symes 
(AKS) theorem, by restricting the ring X(g[ (2)) of Ad* invariants on the dual space qI (2) 
of the full loop algebra to the translate of the coadjoint orbit O//^ C 0[(2)+* C Qi{2)* by the 
fixed element Y. Picking any element (j) G ^aks ^ Hamiltonian, the AKS theorem implies 
that the equations of motion take the Lax form 

^ = [{dcl>){MU,Af] (1.7) 
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and hence determine an isospectral flow. Since ^aks J^^^ ^^^S of spectral invariants, 
fixing simultaneous level sets of its elements amounts to fixing the coeflBcients of the charac- 
teristic polynomial 

nx. := det(A^(A) - a.) = c= - c E + 

i=l ^ ' 



where 



det(A^(A)) -.= ^^ = 1^2 - + M, (1-9) 

n 

a(A) := J](A-ai), (1-10) 



and 



Q / ^ ■ i7 *' 'fc ■ \ 0% 

Z Q!j — OLj \ X\ 

are the (generalized) Devaney-Uhlenbeck invariants (cf. [Mo] ). The leading term of the 
polynomial 

n 

P{X) = Y,P,X' (1.12) 

has constant coefficient 

Pn = -{a' + hc), (1.13) 

while the remaining coefficients {Pq, . . . Pn-i} are independent generators of the ring of com- 
muting invariants. 
Letting 

- y - ^ (1.14a) 

2^A-ai a(A) ^ ^ 

1 ^ 

^:=a(A)(C-^EYzV)' ^1-1^^) 

Z . A QJj 
1=1 



the invariant spectral curve C is defined by 



z'^ = K^{X)-a{X)P{X), 
5 



(1.15) 
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and thus is hyperelliptic. Since -P(A) has leading terms of the form 

/ n " \ 

P(A) = -{a^ + hc)y' + U J] + (aV 6c) J] A^-^ + 0(A"-2), (1.16) 

\ i=l i=l / 

C generically has genus g = n — 1 ii a? + hc ^ ^ or Y^^=i h ^ g = n — 2 ii + be = 
and ELi = 0- 

In the foUowing, instead of considering individual Hamiltonians within the ring of spectral 
invariants, it will be convenient (as in [BT]) to treat the invariant polynomial -P(A) as a 
one parameter linear family of Hamiltonians, all commuting amongst themselves. It will be 
necessary to distinguish two cases, depending on whether c = or c 7^ 0. 

lb. Case (i) c = 0. Restiction to T*5"-^ 

Following the general procedure of [AHH4, H] the spectral Darboux coordinates 
U,pAh,Ci^} 1^=1,. ..n-1 are defined by 



n 



^X-a,- a(A) ^^-^^^^ 

n-1 

Q(A):=n(^-^/^)' (l-17b) 
^x.=l 

n 



and 



e 

i=l 



n 



1 

1=1 

1 

p-=7;^^iVi (1.17e) 



i=l 



Here {g, A^}^=i^...„_i are, essentially, hyperellipsoidal coordinates on R", and {p, C^} ^=i,...n-i 
are the canonically conjugate momenta. To make such an identification, we must also assume 
that the cuj's are all real and positive, and choose an ordering such that, e.g., 

ttn < An-l < Oin-l < Xn-2 < ■ ■ ■ Xl < ai- (1-18) 

The canonical 1-form on R^"^ is 

n n — 1 

= ^ Vidxi = pdq + E Ci^dXij,. (1.19) 

i=l /u=l 
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The Jacobian matrix of the coordinate change is defined by the partial derivatives 

(1.20) 



dXfj, 2 - Q!i ' dq 2 ' 
and its inverse by 

^ 2a(A^) Xj dq ^ 

dxi Q'(A^)A^-aj' dxi 

The coordinate frame fields are thus 

d 1 Xi d 



(1.21) 



1 '^-^ Xi d , , 

= 5ET-37r«;r. /. = i,...n-i (1.22a) 



The R'* euclidean metric in these coordinates is 

n-l 



and the volume form is 



2" inMa(A^)|^ 



dV= ^^'-<^^ ^ ^d^AdAi A...AcZA„_i. (1.24) 



In order to express the linear family of Hamiltonians generated by -P(A) in terms of the 
hyperellipsoidal coordinates, we note first that it follows from (1.6) and (1.8) that the values 
of P(A^) are given by 

a(A^) 4 Xf^-aiJ 



Using Lagrange interpolation, P{X) may then be expressed as 

/ n— 1 n 1 " 

P(A) = Q(A) (-A + J] A^ + ^ «,)«' + ^ E 
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where 

1 h 

-J]/, = 2ap--e^. (1.27) 

i=l 

To write the Hamihon-Jacobi equation, we must reinterpret the coefficients of the invari- 
ant polynomial P{X) in (1.26), not as functions on the phase space, but rather as integration 
constants and replace the canonical momentum components {CfnP}fj,=i,...n-i by the partial 
derivatives {^-, ^}fi=i,...n-i of the Hamilton characteristic function. The resulting form for 
the time independent Hamilton- Jacobi equation is then 

n-l 



Q(A)a(A^) f OS V K{X 



a - 



as 



+ Q(A) ( -a'X + 2a^ - -e^ + a\J2»^-Yl V) ) = (1-28) 



1=1 M=l 

where the leading coefficient Pn of -P(A) is given in (1.13), and the remaining n coefficients 
{Pq, Pn-i} are interpreted as integration constants determining the "energies" for the 

parametric family of Hamiltonians defined in (1.26). 

The integration of (1.28) then proceeds by separation of variables. Expressing S in the 
separated form 

n-l 

S{Xi, . . . Xn-1, q) = So{q) + J2 ^m(Am): (1-29) 

dividing both sides of (1.28) by Q(A) and equating the leading terms in A at oo, as well as the 
residues at {A = A^}^=i^...^_i in the resulting equation gives 

dX^ J a(A^)2 

2a^ = ^e'l -a^Y.ai + P„_i. (1.30b) 

^ i=i 

If a 0, this may be integrated to give the completely separated solution 

n — 1 n 



i=i M=i ^=1*^° V ^ 

(1.31) 

The linearizing coordinates are then 

dS 1 f'^^ A* 
Q.-- =-S^ dX, i = 0,...n-2 (1.32a) 

dP, 2^J, ^K^iX)-a{X)P{X) ' ' 

dS q 1 ^\ f^'^ A"-i 



Qn-1 ■■= ^ = 77- + 7^ V An / dX. (1.32b) 

' dPn-i 2a^2^^ '^io ^K^{X)-a{X)P{X) ^ ' 
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The first n — 1 of these, defined by eq. (1.32a), involve abehan integrals of the first kind and 
essentially define the Abel map to the Jacobi variety J'(C). The last one, defined by (1.32b), is 
an abelian integral of the third kind, the integrand having a pair of simple poles over A = oo. 
The linear flow induced by any Hamiltonian = 0(Po) • • • Pn-i) in the ring I^xs spectral 
invariants is then given by 

Qi = Qio + ^t, i = 0,...n-l. (1.33) 

By standard Jacobi inversion techniques (cf. [Du, GH, AHH4]), any function of the coor- 
dinates {A^, q}^=i,...n-i that is symmetric in the A^'s, can be given an explicit form in terms 
of the Riemann theta functions associated to the curve C. 

If a = and 6 = 0, both q and p are conserved quantities, and Hamilton's equations 
may be integrated on the invariant symplectic submanifold given by fixing a level set of q 
and p. By eqs. (1.17c), (1.1 7e), this defines the cotangent bundle T*S''-'^ C R^" to the 
sphere S'^~^ C R". The Hamilton-Jacobi equation (1.28), may then be interpreted on S^~^, 
by choosing 5" = 5'(Ai, . . . , An-i)- The separated form is again given by (1.29), (1.30a), with 
So = and (1.30b) omitted. Both the leading and next to leading terms in eqs. (1.26) vanish, 
so -P(A) is of degree n — 2. The completely separated solution is 

S(A....,A„_.,) = Ef /^^^^5f^rfA (1.34) 

and the linearizing equations are given by (1.32a) and (1.33), for z = 0, . . . n — 1. Since the 
genus of C is (7 = n — 2, the i = n — 2 integral in (1.32a) becomes singular, the integrand 
having simple poles over A = oo. 

If a = and b ^ 0, q is still a conserved quantity, but p is not. Since the linear family 
of Hamiltonians -P(A) now has no dependence on p, to apply the Hamilton-Jacobi theory, the 

dS 

roles of p and —q must be interchanged, and the term e"^ in (1.28) replaced by e^. The 
Hamilton characteristic function S is now a function of {Xn,p}n=i,...n-i and the solution is 
obtained in completely separated form as 



S(A. .... A„_.p) ^ pin ( + ±t /^-(^)-;W^(^) ,A. (1.35) 



The linear flow equations (1.32a) and (1.33) remain the same for z = 1, ... n — 2, while (1.32b) 
is replaced by 

Q„_i := — = -J— + -\^ dX, i = 0,...n-2. (1.36) 

' aP„_i Pn-i2^J^ ^K^X)-aiX)P{X) ^ ^ 
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Even if a 7^ 0, and q and p are not individually conserved quantities, we may still impose 
the second class constraints 



Y,xl = l, (1.37a) 

n 



or, equivalently, 

g = 0, (1.38a) 
p = 0, (1.38b) 

which define the symplectic submanifold T*S'^~^ C R^" as phase space. The term Pn-i in 
P(A) is no longer viewed as an independent dynamic variable, or an integration constant, but 
rather the fixed constant defined by 

Pn-i :=a2^«i--. (1.39) 

i=l 

The hyperellipsoidal coordinates {A^}^=i,...n-i, given by (1.17a), (1.17b) are now interpreted 
as defined on S'^~^ C R" and, together with the conjugate momenta {Cn} fj.=i....n-i^ these 
provide a canonical system on T*S'^~^. The invariant coefficients {Pq, . . .P„-2} of P{^) still 
form a Poisson commutative set when constrained to T*S'^~^, even though q and p are not 
individually conserved quantities. (This follows from the fact that in the ambient space, the 
combination (1.27) commutes with all the P^'s.) 

For later use in Sec. 2, we note that the induced metric on in the hyperellipsoidal 

coordinates {Ai, . . . A^-i} is 

n-1 



the volume form on S^~^ is 



M — 1 



2" iriu^wi^ 



and the scalar Laplacian is 

n-1 



a(A.) / 1 " 1 d \ , , , 
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The Hamilton-Jacobi equation is the same as in the unconstrained case (1.28), but with 
the ^ term omitted. The solution in completely separated form is again just 



and the linearizing variables and flow are again defined by eqs. (1.32a), (1.33). 

An example of such a constrained integrable system on T*S'^~^ is generated by the 
invariant 

(j)R := -2Pn-2 



n n 

-2a'^aiXiyi + b^aix'^, (1-44) 
i=i i=i 

which, on the constrained manifold defined by (1.37a), (1.37b) becomes 

^R= ^^Vi - -^-^ -"^a^oiiXiyi + b^aiX^. (1.45) 

i=l 3 i=l 1=1 

For a = 0, this gives the Rosochatius system [Mo, GHHW, AHP]. If all the //^'s also vanish, 
it reduces to the Neumann oscillator system [Mo, H]. 

Ic. Case (ii) c ^ Q. Restiction to r*£:"-^ 

In this case, the spectral Darboux coordinates {A^, Cf^}ij,=i,...n are defined by the relations 

V-^ + 2c = 2c^ (1.46a) 



i=l 



n 

Q{X) := HiX - X^) (1.46b) 



n 



1 



The Jacobian matrix is thus again given by 



dX^ 2 Xf^- ai 
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and its inverse by 

dX^ ^ 1 a(A^) Xj 
dxi c Q'{X^) X^- ai' 

The coordinate frame fields are therefore 

d 1 Xi d 



and the canonical 1-form is 

n n 

^ = J]y,da;, = J^CmC^A^- (1-50) 

i=l ij,= l 

To identify {A^}^^i^...„ as hyperellipsoidal coordinates, the constants {ai}i=i^,,,n must again 
be chosen as real and positive, and an ordering fixed, e.g., by 

Xn < an < Xn-1 < OLn-\ < . . . < Ai < Q!i. (1-51) 

The R"^ euclidean metric in these coordinates is 



and the volume form is 

dV = (0 ' A ... A dXn. (1.53) 
For reference in Sec. 2, we note that the scalar Laplacian is 

2^^/^ l^^^y (1.54) 



The values {-P(Aju)}ju=i,...n are again given by 



= - a)^ + ^ (,55) 



and, using Lagrange interpolation, P{X) may be expressed as 



-(A) = - E 1,5^ - - f# ) - «A) (a^ + M . (1.-) 
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To obtain the Hamilton-Jacobi equation, we again reinterpret the coeflBcients of the invari- 
ant polynomial P{X) in (1.56) as integration constants and replace the canonical momentum 
components by the partial derivatives giving 

^ Q(A)a(A,.) (fdS Y K\\S\ „,,w2^,^ p/u n 
" g (A-A„)0'(A,) [[ax;-V - ) - «W (" + M = PW, (1.57) 

where the leading term of -P(A) is — (a^ + 6c) A"' and the remaining terms are independent 
integration constants. 

The integration of (1.57) again proceeds by separation of variables. Expressing S in the 
separated form 

n 

^(Ai,...Aj = J]5^(A^), (1.58) 
dividing both sides of (1.57) by Q(A) and equating the residues at {A = X^}i_i=i,...n gives 

dX^ J a{X^r 
This may be integrated to give 



S(A. .... A„_0 = a E A,. + E / ' ■ ^^W-y^W rfA. (1.60) 

The linearizing coordinates are thus 

Qi:=^ = -Y dX, i = 0,...n-l. (1.61) 

dPi 2^iio V^2(A)-a(A)P(A) ^ ^ 

The first n — 1 abelian integrals in (1.61), with i = 0, . . .n — 2, are all of the first kind, again 
defining the Abel map to the Jacobi variety J{C), while the remaining one, giving Qn-i, 
is again singular, the integrand having a pair of simple poles over A = oo. The linear flow 
induced by any Hamiltonian (f) = 4){Po, . . . Pn-i) in the ring of spectral invariants I^ks ^ 
before, given by 

Q^ = <5^o + |^^, z = 0,...n-l. (1.62) 
We may also impose the constraints 

n 2 

J2— = ^c (1.63a) 

n 

J2^ = (1.63b) 



1=1 
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defining the cotangent bundle T*£^"~^ to the eUipsoid £'^~-^ C R"' defined by (1.63a). (Here, 
the constant c must be chosen as positive; more generally, arbitrary signs may be allowed for 
the Qfi's and c, thereby defining various hyperboloids.) In terms of the canonical coordinates 
Cij,}ix=i,...n, the two constraints are equivalent to 



An=0 
Cn = 0. 



(1.64a) 
(1.64b) 



For reference in Sec. 2, we note that the induced metric on S"' is 



i=l 



(1.65) 



the corresponding volume form is 



-1 n— 1 



n-1 



£n-l 



(1.66) 



and the scalar Laplacian is 



2 



1 

+ 



n 

— - 



_d_ 

dX, 



(1.67) 



Eqs. (1.64a), (1.64b) determine the cotangent bundle T*S'^~^ C R^" as a symplectic 
submanifold; i.e., they are purely second class, and the restriction of the remaining coordinates 
{AjU, Cij,}ij.=i,...n-i provide canonical coordinates on this constrained manifold. The restriction 
of the canonical 1-form defines the canonical 1-form on T*S^~^ 



n-l 



9\T*£ri-i — C^dXij,. 



(1.68) 



Although neither of the constraints (1.63a), (1.63b) is individually invariant under the AKS 
fiows generated by the invariants Pq, . . . Pn-i, they are equivalent to the pair 



Q(0) = 
a(0) a(0) ^ 4 [^^ 



Oli 



(1.69a) 
(1.69b) 
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and the second of these is invariant. It follows that the restrictions of the remaining invariants 
{Pi, . . . Pn-i} to the constrained manifold T*E'^~^ C R^" also Poisson commute, generating 
completely integrable systems. On T*S'^~^, we may write 

g(A) = AQo(A) (1.70) 

where 

n-l 

QoiX)--=l[{X-X^). (1-71) 
Proceeding again by Lagrange interpolation, we have 

+ go(A)(-(a^ + .c)A + ^(-a^ + i(|:g)^]], (1.72) 



which, with suitable reinterpretation of P{X) in terms of integration constants, gives the 
Hamilton-Jacobi equation on T*£'^~^ as 

""'^ Aa(A^) ( ( dS V K^{X, 



+ 0„(A) {-(a- + MA + ^ i-a- + i (t £ ] ) ) = PW, (1.73) 



with P(0) fixed to satisfy (1.69b). Expressing S'(Ai, . . . A„_i) in separated form as 

n-l 

5(Ai,...A„_i) = ^^^(A^), (1.74) 

and proceeding as in the unconstrained case again gives (1.59), for // = l,...n — 1, which 
upon integration gives 



5(A....,A„_.) = a"fA, + g f' I'm^J^i,. (l.rs) 



The linearizing coordinates are therefore 



Q..- =-\^ dX, i = l,...n-l. (1.76) 
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and the linear flow induced by any Hamiltonian = 0(Pi, . . . P„_i) in the ring T\j^g of 
spectral invariants is again given by 

Qi = Qio + ^t, z = l,...n-l. (1.77) 

In general, the abelian integrals in (1.76) are all of the first kind, so the linearizing map is 
again the Abel map to the Jacobi variety J^{C). 

As an example of such a constrained system within the commuting family of systems 
generated by ^aksi consider the Hamiltonian 



(f)j := 2Pn-i = ^ + (a^ + be) «i 

i=l ^ « ^ i=l i=l 

Taking c = ^, a = 6 = 0, {jii = 0}j=i,...n, this reduces to the Jacobi system, determining 
geodesic motion on an ellipsoid [Mo]. More generally, if 6 7^ 0, /ij 7^ 0, we have motion on 

the ellipsoid in the presence of both a harmonic force, directed towards the origin, and forces 

2 

derived from the Ym=i ^ potential, as in the Rosochatius system. 

2. Quantum Integrable Systems in gl(2)+* 

2a. Ambient Space. 

We quantize the systems considered above first in the ambient space R", using canonical 
quantization within the Schrodinger representation. The wave functions will thus be taken as 
smooth functions ^^(0:1, . . . Xn) on R"^, square integrable with respect to the standard measure. 
Using the canonical quantization rule (with ^ = 1) 

d 

Xj — >x:= Xj Vj — > Vj := -i— := -idj, (2.1) 
we obtain a quantum analogue of the loop algebra element Af{X) 

where 

h{X) Yl \_ + 2a (2.3a) 
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^2 , M,^ 



^(^) E T377^ + 26 (2.3b) 
j=i •? 



n 2 
X 



/^(^) ^= E + 2c (2.3c) 



^(^) ^= E -13^ - 2«- (2-3d) 

Note that we have not chosen to order these operators so as to assure self-adjointness; they 
wiU only enter at an intermediate stage in subsequent calculations and need not themselves 
be viewed as quantized dynamic variables. We also define an operator that plays essentially 
the role of detA/'(A) in the preceding calculations: 

P(A):=/i(A)^(A)-/(A)e(A). (2.4) 

Again, we are not concerned with whether this is an appropriately ordered "quantum deter- 
minant", since it will only appear as a calculational convenience in what follows; the choice 
of ordering in (2.4) is the most convenient in the subsequent Lagrange interpolation. 
A direct computation shows that I^(A) may be expressed as 

^ n f ^ n n ^ l"- Q . n ^ 

(2.6) 

where 



Ii:=^Y. ' ' 2mOr,a, + - -K'-c + ^ 

2^ ai-aj 2 V 

(2.6) 

Here the /^'s are self-adjoint, and represent the quantized version of the Poisson commuting 
invariants li of the preceding section. Moreover, a direct computation shows that these /j's 
also commute amongst themselves 

[7i,/,-] = V z,j = l,...n. (2.7) 

Using either definition (1.17a)- (1.17b) or (1.46a)-(1.46b) for the coordinates A^, and the 
fact that /(A^) = 0, it follows from eq. (1.22a) or (1.49) that, 

^(A.) - ^ - 2za^ + - ^ ^^--^ + - a , (2.8) 
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where it is understood that the factor (A^ ~ in the denominator precedes the derivation 
Xjdj. Define the difi^erential operator valued polynomial P(A) by 



a(A) ■ 2 ^ A - a 



n f 



Equating the expressions for T>{\ij) given by eqs. (2.5) and (2.8), gives 

a(A^) V^A^ / 2^A/,-a, V^A^ / o2(A/,) 

where again, the ordering on the LHS of (2.10) is understood to place the a(A^) term in the 
denominator before the difi^erential operator in the numerator. Note also, that in the definition 
(2.9) of -P(A), the leading term is 

1 

P{\) = ^Y.^Iiy'-^+0{\n-2\ (2.11) 



2 



where 



4-1 := 2 E = - 2 E 2^a(^»5^ + 2 ) + + ^(^' + §) ) • (2-12) 

Thus, contrary to the definition of -P(A) in (1.9), -P(A) here does not contain the constant 
leading term — (a^ + bc)X^, and is an operator- valued polynomial of degree n — 1. 

In the following sections, we again must distinguish between the cases c — and c 0. 

2b. Case (i) c = 0. Restiction to 5^"^ 

Taking (1.17a)-(1.17c) as the definition of our coordinates {A^, q}ij,=i,...,n-i ^ind of Q(A), 
and using eqs. (2.9), (2.11), (2.12), Lagrange interpolation gives 

A.xN Q(A)a(AJ (fd .A'lA 1 fd .\ K\X.) 



as the quantum analogue of (1.26). Here we have used 
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Letting 

^-ItJL- (2 16) 



where {/f'}i=i,...n is a set of simultaneous eigenvalues of the operators {/i}i=i,...n, the time in- 
dependent Schrodinger equation corresponding to the 1-parameter linear family of commuting 
Hamiltonians -P(A) is thus 



//o \2 



d nia b 
''dq~ 2 



+ Q{X) ( -2ia§- - ^ - I E * = P^(A)*, (2.16) 



where ^ = ^{Xi, . . . X^-i, q). 

Choosing \E' in the completely separated form, 



n-l 



- v-i, n-l, - rw; 11 rMV-M/- ''=^ ^ (2-17) 



dividing (2.16) by Q(A) and equating the leading terms in A and the residues at A = A^, we 
find that each V'm('^m) satisfies the same separated equation 



while (f){q) satisfies 



(-2za| - ^ - ^e'^) </>(,) = i f: /f 0(,). (2.18b) 



Eq. (2.18a) is a Fuchsian differential equation of the generalized Lame type already encoun- 
tered in previously studied examples of such systems [BT, Mc]. Note also that, if A and ( are 
viewed as canonically conjugate variables, (2.18a) may be interpreted as the quantized form of 
the characteristic equation (1.15) defining the classical spectral curve C. Thus, the completely 
separated form of the Schrodinger equation may be viewed as the "quantized spectral curve" 
(cf. [Sk3]). 

If a = 0, g is a conserved quantity and if 6 7^ 0, eq. (2.18b) should be interpreted as fixing 
the constant value of q in terms of the invariant ^"^^ if ■ Then eq. (2.18b) implies that, 
rather than choosing (p to be an function, it should be chosen as a delta function within 
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the coordinate representation in which q is fixed by the eigenvalues if . If a = and 6 = 0, 



this choice may still be made, but the eigenvalues if are not all independent, since they must 



sum to 0, 

n 

Y^lf^^, (2.19) 

but the invariant q may be set equal to any constant value independently. 

Whether a = or not, we may impose the second class constraints (1.38a), (1.38b) in the 
quantum problem by requiring the wave function to satisfy 

gj + -j*|,.„ = 0, (2.20) 

which is the quantum analogue of the constraint (1.38b), and setting = in the expression 
(2.13) defining -P(A). Restricting to the subspace consisting of ^^'s satisfying (2.20), the 
coefficients of P{X) generate an n — 1 dimension linear space of commuting operators on the 
sphere S'^~^. This is spanned, e.g., by the restrictions of the operators which, because of 
(2.20) satisfy 

n 

Y.ij^-h. (2.21) 
Applying the constraint, (2.16) reduces to 

-6g(A)*o = i'^(A)*o, (2.22) 

where 

nq 

*(Ai,...An-i,g) :=e--*o(Ai,...An-i) (2.23) 
and the eigenvalues if must satisfy 

n 

Y.lf = -h. (2.24) 

The volume form on S'^~^ in these coordinates is given by (1.41) and the polynomial family of 
operators appearing on the LHS of (2.22) are all self-adjoint with respect to the corresponding 
scalar product. The separation of variables for the resulting systems on the constrained space 
S'^~^ is again obtained by setting 

n— 1 . . , 

*o(Ai,...A,_i) = nV;^(A^)e^"^M=i^'^ (2.25) 

i=l 
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in (2.22), dividing by Q{X) and equating residues at A = A^. The resulting separated equations 
for the ipi_i{Xfj,) are all again given by (2.18a). The particular case of this family of systems 
corresponding to choosing all //j = and a = 0, as in the classical case, gives the quantum 
Neumann oscillator system (cf. [BT]), while if the //j's are non-zero, we obtain the quantum 
Rosochatius system (cf. [Mc]). 

2c. Case (ii) c^O. Restiction to £'^~^. 

We begin again with the unconstrained system in the ambient space R". Taking (1.46a), 
(1.46b) as the definition of the coordinates {A^}^=i^...„, and of (5(A), and again using eqs. 
(2.9)-(2.12), Lagrange interpolation gives 

P(A) = (2.26) 



^ Q(A)a(A^) l(d . V ,1 1 (d .\^K\\^) 

2^ (\-\. \r)'(\. \\\ ^o2^\..- 



— ia \ -\ „,\ ^/ + + bc 



(A - A^)g'(A^) y\dX^ J 2 A^ - \dX^ J a^{X,) 

as the quantum analogue of (1.56). Defining the polynomial family of eigenvalues P^{X) 
again by eq. (2.15) in terms of the set {/j^ji^i ...„of simultaneous eigenvalues of the operators 
{Ii}i=i,...n, the time independent Schrodinger equation corresponding to the l~parameter 
linear family of commuting Hamiltonians -P(A) is now 

Q(A)a(AJ [fd .Yl^ 1 fd .\ K^iX.) 2 
(A-A^)Q'(A^) \^V^A^ / 2 A^ - V / «'(Am) J 

= P^(A)*, (2.27) 
where = ^'(Ai, . . . A^)- Choosing in the completely separated form, 

n 

*(Ai, . . . A„) = e'''^l=^ H V'^(A^), (2.28) 

dividing (2.16) by Q{X) and equating the residues at A = A^, we find again that each il^n{X^) 
satisfies the same separated equation 

.7 = 1 ■' ^ ' .7 = 1 



In the the quantum problem corresponding to the system with second class constraints 
(1.64a) (1.64b), it is preferable to impose the invariant constraint (cf. eq. (1.69b)) 



P(O)^' = P^(0)*, (2.30) 
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where 



Letting 



a(0) ~~2^ ai ~ a2(0) 



+ bc. 



_d_ 
dX, 



la 



1^ 1 



d 



ia + 



+ + he, 



we have 



n-l 



(A-An)Qo(A)a(A^) 



Qo{Xn) 

Qo(A) 



where 



^(A) - ^(A.) + A, 

a(A^) 



n-l 



B(A):=E7T3 



;^(A-A„Wi(AJ " 
It follows that, for satisfying the constraint (2.30), we have 



P(A)* = Po(A)*, 



where 



PoW := Qo{X) R{X) - R{0) + 



P^(Q) 
Qo(0) 



= Qo(A) 



n-l 



Xa{Xf,] 



^^{x-x^)XMM 



P^(0) 

Qo(o) 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



Since the coordinate An does not enter in eq. (2.35), it may be treated as a parameter, 
and set equal to in This may be viewed as the gauge condition associated to the first 
class constraint (2.30). Since the constraint commutes with the family of operators -P(A), 
the reduced operators -Po(A) defined by (2.36) still commute amongst themselves. The wave 
function with A^ = 0, is defined on the ellipsoid S^~'^, and satisfies the polynomial family 
of Schrodinger equations 



go(A) 



n-l 



Aa(A^) 



P^(0) 



]^^{X-X,)XMM^^ ' Qo(0) 



* = P^(A)*, 



(2.37) 
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where P^(0) is given by (2.31). The coefficients of -Po(A) thus generate an n — 1-parameter 
family of commuting invariants, which may be simultaneously diagonalized. The simultaneous 
eigenfunction is independent of the parameter A, so dividing (2.37) by Qo{X) and taking 
the limit A — > oo gives 

and hence 

Po{X) = Qo(A) (^R{X) + J2 j (2-39) 
The Schrodinger equation (2.37) may therefore equivalently be written 



go(A) 



n— 1 / \ \ n ' 

^^(A-A,)A,Q^(A,) 



^ = P^(A)*. (2.40) 



Either way, ^'(Ai, . . . A^_i) may be chosen in completely separated form 



n-l 



Ert — 1 

^y..„.....n-LJ-^ " n^M(^M)> (2-41) 

where each ijj^{Xfj_) still satisfies (2.29). 

Notice, however, that the reduced operator Pq{X) obtained in this way is not self-adjoint 
with respect to the scalar product 

:= / §^dVsr.-i, (2.42) 

with dVsn-i the standard volume form (1.66) on £'^~^. It is, however, with respect to the 
scalar product 

(^7*) / (2.43) 



where dV£n-i is the modified volume form defined by 

n—i n—1 n—1 



Making the transformation 

n-l 
/u=l 

n—1 n—1 

Po ^ ^0 := n ° ^0 ° n (2.45) 

11=1 /i=i 
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gives the equivalent Schrodinger equation 



Po(A)* = P^(A)*, 
where the polynomial family of operators Po{X), given by 



(2.46) 



n-l 



Po(A):=E7T^ 



Ag(A)a(A, 



_d_ 



la 



1 J_ 



+ 



«'(Am) 



+ + bc + U{X^) 



_d_ 
dX, 



— la 



(2.47) 



with 



do 



(2.48) 



is now self-adjoint with respect to (2.42). 

The terms IA{X^) in (2.47) are of order 0{fi^), and hence disappear in the semi-classical 
limit. If we subtract these, we obtain the polynomial family of self adjoint operators 



n-l 



^c(A) := J2 



XQ{X)aiX, 



(A - A^)A^g'(A^) 



_d_ 
dX, 



ia + 



^■^1 A^ - aj 



A, 



_d_ 
dX, 



la 



+ 



K\X,) 
«^(Am) 



(2.49) 



which have the same semi-classical limit as (2.47). The coefficients of -Pc(A) may again be 
simultaneously diagonalized by separation of variables in the corresponding Schrodinger equa- 
tion 

Pe(A)*c = i'c^(A)*c, (2.50) 



where the eigenvalues defined by 

^c^(A) 
a(A) 

are again constrained to satisfy 



jE,c 



A — CKi ' 



(2.51) 



Pf(0) 



n jE,c 



9 , 



K2(0) 



a(0) 2^ tti a2(0) 
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Expressing in the completely separated form 

E, n— 1 
71 — 1 , 

"=^"11^-.^^, (2-53) 

each il^cixiXfj) must satisfy 

Since the self-adjoint operators Po{\) and -Pc(A) have the same semi-classical limits, and 
both are completely integrable through separation of variables, either may be regarded as a 
valid quantization of the corresponding constrained classical system on £^~^ . Ambiguities of 
order OiTi^) are well known to occur in the constrained quantization procedure [D], and can 
only be resolved by appealing to some further physical principle. In the case a = 6 = //j = 0, 
the degree n — l coefficient in -Pc(A), which in the classical system generates free motion on 
£'^~^ (cf. eq. (1.78)), gives the Laplacian on £^"~-^, whereas that in Pq{\) contains additional 
Oijhb'^) potential terms due to the W(A^) factors in (2.47). This may be a valid reason to regard 
Pc(A) as the correct quantization of this family of integrable systems. 

3. Discussion 

We have seen that in the ambient space R"^ and on the sphere, the quantization of the in- 
tegrable systems of Section 1 proceeds straightforwardly, whereas on the ellipsoid, constrained 
quantization leads to 0{Tt?) ambiguities. The origin of the problem may be understood by 
noting that, whereas in the coordinates of Sections lb, 2b, which are adapted to the sphere, 
not only does separation with respect to the transversal coordinate q hold for the ambient 
space Schrodinger equation, but also for the volume form (1.24), giving rise to the standard 
sphere volume form (1.41) upon factorization. In the case of the ellipsoid, however, although 
the Schrodinger equation again separates in the adapted coordinates of Sections Ic, 2c, the 
volume form (1.53) does not admit such a factorization with respect to the transversal coor- 
dinate A„. If one just imposes the constraint A„ = defining the ellipsoid, and eliminates 
the dXn factor in (1.53), the volume form (iV^n-i of eq. (2.44) results, not the standard one 
(1.66) determined by the induced metric. More generally, when the constrained manifold is 
obtained as an orbit in configuration space under a large symmetry group, and the commut- 
ing operators of the system are expressible, as for the sphere, in terms of invariant operators 
in the enveloping algebra, the constraining procedure may compatibly be deduced through 
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separation of the remaining "radial" variable. This is the case, for example, in the case of free 
motion on symmetric spaces [BKWl, BKW2, K, KM, ORW]. For the ellipsoid, or more 
general constrained systems in which no such transitive symmetry group is present, further 
considerations regarding self-adjointness invariably enter, and 0{?i^) ambiguities of the type 
encountered here may arise. 

The systems studied here (prior to the constraints) may, in view of the form of the op- 
erators A/'(A) of eq. (2.2), and the choice of commuting invariants given by its "quantum 
determinant" , be regarded equivalently as defining a noncompact version of the Gaudin spin 
lattice [G, Ski— Sk3, Ku] (with su{2) replaced by s/(2, R)). From this viewpoint, the sepa- 
ration of variables may be interpreted as a "functional Bethe ansatz" (cf. [Ski— Sk3]). 

Two types of generalization of the systems studied here naturally suggest themselves, 
both at the classical and quantum levels. The first is to allow degeneration of the parameters 
{ai, . . .an} defining the image of the Poisson map (1.1b), (1.2), including the possibility 
of higher order poles. This is known (cf. [K, Ku]) to lead to separable coordinates other 
than the generic, hyperellipsoidal ones appearing here, adapted to symmetry groups larger 
than the finite group (Z2)"' encountered in Section 1. The second generalization consists of 
extending the present considerations to higher rank algebras, such as £|((r), and reductions 
thereof. Whereas the classical systems so arising are known to be integrable and separable 
under generic assumptions regarding the rational coadjoint orbits and initial data [AHH4], 
very little is known about the corresponding quantum systems. These questions merit further 
study, and will be addressed in future work. 

Acknowledgements. The authors have benefited from helpful discussions with C. Duval, E. 
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